N9 2- 15 864 


1991 

NASA/ASEE SUMMER FACULTY FELLOWSHIP PROGRAM 

MARSHALL SPACE FLIGHT CENTER 
THE UNIVERSITY OF ALABAMA IN HUNTSVILLE 


IGNITION TRANSIENT ANALYSIS 

Prepared By: 

Academic Rank: 

Institution: 

NASA/MSFC : 

Laboratory: 

Division: 

Branch: 

MSFC Colleague: 

Contract No . : 


OF SOLID ROCKET MOTOR 


Samuel S. Han, Ph.D. 

Associate Professor 

Tennessee Technological 
University, Department of 
Mechanical Engineering 


Propulsion 
Propulsion Systems 
Performance Analysis 

S. Don Bai, Ph.D. 

Charles F. Schafer, Ph.D. 

NGT-01-008-021 

The University of Alabama 

in Huntsville 




INTRODUCTION 

Measurement data on the performance of space shuttle solid 
rocket motors show wide variations in the head-end pressure 
changes and the total thrust build-up during the ignition transient 
periods. To analyse the flow and thermal behavior in the tested 
solid rocket motors, a 1-dimensional, ideal gas flow model via 
SIMPLE algorithm was developed [4]. Numerical results showed that 
burning patterns in the star-shaped head-end segment of the 
propellent and the erosive burning rate are two important factors 
controlling the ignition transients. 

The objective of the present study is to extend the model to 
include the effects of aluminum particle commonly used in solid 
propellents. To treat the effects of aluminum-oxide particles in 
the combustion gas, conservation of mass, momentum and energy 
equations for the particles are added in the numerical formulation 
and integrated by following IPSA [5] approach. 


METHOD OF ANALYSIS 


Governing Equations 

Conservation of mass is 
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where subscript i indicates i component, A is the cross-sectional 
area of the port, r w is the volume fraction , P is the material 
density, u is the velocity, is the mass fraction in the 
propellent, is the density of the propellent, b is the perimeter 
of the burning cross-section, is the igniter gas flow rate. 
Burning rate of the propellent is given by [2] 


= r rt/ c f/r ief ) . (2) 

The first part in the RHS of eqn (2) is the standard burning rate 
and the second part the erosive burning rate. Conservation of 
mass also requires that 
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Conservation of linear momentum is 
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Thermodynamic pressure of the system is due to gas phase only 
and is given by 


P = A T > / IP, ) , (5) 

where f ( P, ) is a correction needed for high pressure gas [3], The 
wall shear stress is given by 
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where f is the friction factor obtained by measurement. 
The inter-phase friction force is given by [1] 
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where d % is the particle diameter and the Reynolds number is 
defined by 
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Energy conservation is , in terms of internal energy. 
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where k’ is thermal conductivity, h^ is the enthalpy of formation, 
h;, is the enthalpy of igniter gas, and ±i is the dissipation. 

Heat transfer from the gas to the solid propellent when the 
propellent surface temperature is lower than a preset auto-ignition 
temperature is determined by 

l - K t T - T, (10) 

where h c is a convective heat transfer coefficient determined by 
measurement. Inter-phase heat transfer is calculated by [1] 
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where 
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Heat transfer from the gas to the solid propellent is assumed 
to be one dimensional transient conduction [4]. 

Initial and Boundary Conditions 

Initially stagnant atmospheric air is filled in the chamber 
and the nozzle. Transient process begins with the introduction of 
igniter gas into the chamber at the head— end section. Heat transfer 
from the gas to the solid propellent raises temperature of the 
propellent resulting in auto-ignition. Mass and energy released 
from the burning solid propellent rapidly increase momentum and 
energy of the flow in the system and the total thrust. Open 
boundary conditions are maintained at the exit of the nozzle and 
solid wall conditions are used at the head-end. 


N-2 



Numerical Method 

IPSA (Inter-Phase-Slip-Algorithm) originated by Spalding [5] 
is used to approximate the solution of governing equations. IPSA is 
an extension of SIMPLE method and consequently follows a similar 
computational procedures. Equation (1) is solved for the volume 
fraction of particle phase and equation (3) is used to find volume 
fraction of the gas phase. With known volume fractions and a 
guessed pressure field, momentum equations are solved for 
velocities. Corrections on density, velocity and pressure are made 
by solving pressure correction equation. Energy equation is solved 
to update temperature change. These steps are repeated until 
convergence is satisfied within a given time step. 

RESULTS AND CONCLUSIONS 

Geometry of the motor, igniter gas flow rate and other 
physical and numerical parameters are those used in a previous 
analysis based on a single fluid model [4], 

Numerical results obtaind by the single fluid option of the 
present model are shown in Figure 1 in comparison with measurement 
data. Numerically predicted pressure increase after the flame 
spreding (0.12 sec) is higher than the measured data. Higher 
pressure increase is due to higher burning rate which is caused by 
higher flow speed in the chamber. Predicted total thrust is much 
higher than the measured value (about 46 % higher at 0.6 sec). 

Numerical results using two-fluid option of the model are 
shown in Figure 2 . Particle diameter is assumed to be 10 m and the 
density of the particle is 1500 kg/m . Particle loading of the 
aluminum in the solid propellent is 10 %. Due to increased mixture 
density, flow speed in the chamber is much slower and subsequently 
slower increases in pressure gradient and total thrust. Numerical 
results are in agreement with physical expectations. 
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